
function zout = zshape(r,theta)

% Elliptical domains (Huisinga).

% Shape of the domain. The Leja points will be distributed over the 
% of the domain "bondary". 
 % This shape correspond to an ellipse centered -r.
      zout =cos(theta)*(1+r)-(1+r)+i*sin(theta)*(1.d0-r); 
      end
